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Ground state and thermodynamics of geometrically frustrated spin-1/2 Ising-Heisenberg model 
on two different but topologically related triangles-in-triangles lattices is investigated in particular. 
A rigorous mapping based on generalized star-triangle transformations enables one to derive exact 
analytical results for the free and internal energy, spontaneous magnetization, entropy and specific 
heat from the known exact results for the corresponding spin-1/2 Ising model on a triangular lat- 
tice. It is demonstrated that several unconventional quantum ordered and disordered phases may 
become the ground state in addition to the classical ferromagnetic phase with a perfect alignment 
of all constituent spins. A mutual competition between the ferromagnetic Ising and XXZ Heisen- 
berg interactions may lead to an appearance of a peculiar quantum ferromagnetic phase, in which 
the magnetic moment of the Heisenberg spins undergoes a quantum reduction to one third of the 
saturation value in contrast with the fully saturated magnetic moment of the classical Ising spins. 
In the highly frustrated regime, we have found either the disordered quantum paramagnetic phase 
with the substantial residual entropy or a similar but strikingly spontaneously long-range ordered 
phase, in which a quantum reduction of the magnetic moment of the Heisenberg spins causes a small 
quantum reduction of the magnetic moment of otherwise purely classical Ising spins. It is evidenced 
that the specific heat exhibits diverse thermal dependences with or without a logarithmic divergence 
superimposed either on one or two round Schottky-type maxima. 
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I. INTRODUCTION 



Frustrated two-dimensional (2D) quantum spin mod- 
els belong to the most challenging issues to deal with in 
the modern statistical and condensed matter physics^— 
A considerable interest in investigating 2D geometri- 
cally frustrated spin models is closely connected with 
an existence of numerous layered insulating magnetic 
materials, 4-7 which can be characterized by incapabil- 
ity of spins to simultaneously satisfy all pair spin-spin 
interactions. In general, the magnetism of the vast ma- 
jority of insulating materials is well captured by the 
quantum Heisenberg model and its various extensions. 8 
It should be mentioned, however, that a mutual inter- 
play between the geometric frustration and quantum 
fluctuations is not only responsible for an extraordi- 
nary diverse magnetic behavior, but is also the main 
cause of a mathematical intractability of geometrically 
frustrated Heisenberg models with exception of a few 
prominent exactly solved cases^ On the other hand, 
the classical Ising model is still exactly solvable on sev- 
eral frustrated 2D lattices and it may thus provide a 
valuable insight into cooperative phenomena originating 
from the spin-frustration effect j 10 ' 11 The most limiting 
drawback of the Ising model unfortunately lies in the 
lack of appropriate experimental realizations which are 
rather scarce^^ Among the most striking phenomena 
emerging in geometrically frustrated classical and quan- 
tum spin systems one could mention the existence of 



reentrant phase transitions ; 10 ' 11 the quantum reduction 
of the magnetization; 13 ' 14 the presence of magnetization 
jumps and intermediate plateaus in the magnetization 
process ; 3 ' 15 : 16 " order- from-disorder" effect ; 17 : 18 the ap- 
pearance of several unusual spin-liquid ground statesji^^ 
the enhanced magnetocaloric effect ji2r—, the localized 
magnons in a close vicinity of saturation fielder— etc. 

Despite a great effort in this research area, the classi- 
cal and quantum frustrated spin models defined on 2D 
" triangles-in-triangles" (TIT) lattices have received much 
less attention so far. TIT lattices generally consists of 
smaller triangular entities embedded in either some or 
all triangular cells of 2D triangle-based lattices such as 
triangular or kagome lattice. The sites of original lat- 
tice might be thus called as nodal lattice sites, while 
the sites of smaller triangular entities may be viewed as 
decorating lattice sites. In this way, one may consider 
various geometrically frustrated spin models on differ- 
ent but geometrically related TIT lattices. If all lattice 
sites are occupied by the Ising spins, one considers the 
classical Ising model whose magnetic behavior is not af- 
fected by quantum fluctuations at all. If all lattice sites 
are occupied by quantum Heisenberg spins, one turns to 
the quantum Heisenberg model that is basically affected 
through quantum effects. If the nodal lattice sites are 
occupied by the Ising spins and the decorating sites by 
the quantum Heisenberg ones, one considers the hybrid 
Ising-Heisenberg model affected through local quantum 
fluctuations only. 

Up to now, there exist only a few theoretical studies 
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of classical and quantum spin models on one very special 
example of TIT lattice called as the triangulated kagome 
lattice^ - — A particular interest in studying the role of 
geometric frustration in the triangulated kagome lattice 
has been stimulated by the family of polymeric coordi- 
nation compounds CugX2(cpa)6-nH20 (X=F,Cl,Br and 
cpa=carboxypentonic acid), which provides an experi- 
mental realization of highly frustrated magnetic materi- 
als with a remarkably high frustration ratio (> 130) and 
possibly spin-liquid ground state.— ~— In addition, there 
are also experimental indications of one-third plateau 
in the low-temperature magnetization process ) 40 ' 41 while 
the temperature dependence of the susceptibility implies 
two basically different exchange interactions in this in- 
triguing series of magnetic compounds.— ~— A stronger 
antiferromagnetic coupling has been assigned to the ex- 
change interaction inside the smaller triangles and a 
possibly much weaker ferromagnetic coupling to the ex- 
change interaction within the larger triangles of the tri- 
angulated kagome lattice i 40 ' 41 ESR measurements gave 
the almost isotropic g-value and this result is taken to 
mean that CugX2(cpa)g-nH2 compounds afford a suit- 
able experimental realization of the isotropic Heisenberg 
model i 39 ' 41 Of course, the quantum Heisenberg model 
on the triangulated kagome lattice cannot be rigorously 
solved and thus, this model has been dealt with the help 
of different approximate perturbation approache a 31 ' 32 
and numerical exact diagonalization^ 3 - - — However, the 
substantial difference between the exchange interactions 
within the smaller and larger triangles of the triangu- 
lated kagome lattice provides a plausible ground for intro- 
ducing the simplified Ising-Heisenberg model, which cor- 
rectly reproduces the stronger Heisenberg coupling inside 
the smaller triangles and merely approximates the much 
weaker exchange interaction inside the larger triangles by 
the Ising coupling. A comparison between the exact so- 
lutions of the Ising 26 ' 27 and Ising-Heisenberg 2 ^— models 
on triangulated kagome lattice reveals that a presence of 
local quantum fluctuations leads to a marked decrease in 
the residual entropy of the quantum paramagnetic state 
of the latter classical-quantum model. 

On the other hand, there is a lack of theoretical studies 
concerning with the classical and quantum spin models 
defined on other geometrically frustrated TIT lattices, 
which could be derived for instance from a simple trian- 
gular lattice. The main goal of the present work is to 
fill in this gap, whereas the main emphasis will be laid 
on a comparative study of the Ising-Heisenberg models 
defined on two different but geometrically related TIT 
lattices each of them descended from a simple triangu- 
lar lattice. It will be demonstrated that two different 
magnetic architectures, even though they are geometri- 
cally closely related, will cause a different level of local 
quantum fluctuations and thus fundamentally change the 
magnetic behavior of the considered spin models espe- 
cially in the highly frustrated regime. It actually turns 
out that the ground state of the first considered model 
in the highly frustrated regime is disordered in contrast 



to the partially ordered and partially disordered ground 
state of the second model, which is characterized by the 
obvious quantum reduction of the spontaneous magneti- 
zation. 

The organization of this paper is as follows. Both in- 
vestigated models are introduced in Section |TT] along with 
basic steps of their exact analytical treatment. The most 
interesting results for the ground state, order parameters 
and overall thermodynamics are then discussed in Sec- 
tion IIIII Finally, several concluding remarks are drawn 
in Section HVl 



II. MODEL AND SOLUTION 

Let us introduce the spin-1/2 Ising-Heisenberg model 
on two different TIT lattices, which are schematically 
depicted on the left-hand-side of Fig. QTa)-(b). As il- 
lustrated, the Ising spins a — 1/2 placed at nodal lat- 
tice sites of a simple triangular lattice (full circles) are 
mutually inter-connected through smaller triangles of 
the quantum Heisenberg spins S = 1/2 (empty circles), 
which are either placed into each up-pointing triangle 
(Fig-da)) or into each triangle (Fig.rTfb)) of the underly- 
ing triangular lattice. From this perspective, the smaller 
triangles of three quantum Heisenberg spins (Heisenberg 
trimers) are embedded in the larger triangles of the clas- 
sical Ising spins and hence, both investigated magnetic 
structures fall into the class of TIT lattices. The total 
Hamiltonian of the spin-1/2 Ising-Heisenberg model on 
both aforedescribed TIT lattices can be defined as 
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where a%, and Sf (a — x, y, z) denote spatial compo- 
nents of the spin-1/2 operator of the Ising and Heisen- 
berg spins, respectively. The parameter Jh stands for 
the XXZ Heisenberg interaction between the nearest- 
neighbor Heisenberg spins from the small inner triangles, 
A is a spatial anisotropy in this interaction and the pa- 
rameter J\ represents the Ising-type interaction between 
the nearest-neighbor Ising and Heisenberg spins. 

For further convenience, let us rewrite the total Hamil- 
tonian as a sum over the six-spin cluster Hamiltonians 
schematically illustrated in Fig. [Uc)-(d) 



u 



k=l 



(2) 



Here, N stands for the total number of the Ising spins 
and 7^ denotes the total number of six-spin clusters, i.e. 
7=1 for the TIT lattice drawn in Fig. QJa) and 7 = 2 
for the TIT lattice shown in Fig. [IJb). It is of principal 
importance that all the interaction terms including the 
three Heisenberg spins from the fcth inner triangle belong 
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FIG. 1: (a)-(b) The spin-1/2 Ising-Heisenberg model on two different but geometrically related TIT lattices and its rigorous 
mapping to the spin-1/2 Ising model on a triangular lattice. The full circles denote lattice positions of the Ising spins a = 1/2 
and the empty ones lattice positions of the Heisenberg spins S — 1/2; (c)-(d) Schematic representation of two equivalent 
star-triangle transformations used for establishing a precise mapping correspondence between both the models. 



to the same cluster Hamiltonian %k 
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which in turn ensures a validity of the standard com- 
mutation relation [Hi,'Hj] — between different cluster 
Hamiltonians i ^ j. Owing to this fact, the total parti- 
tion function can be partially factorized into a product 
of cluster partition functions and one may trace out spin 
degrees of freedom of different Heisenberg trimers inde- 
pendently of each other according to the relation 

■yN 7 AT 

z = e n Trfe c M-pH k ) =j2ii Zk ( 4 ) 

{<Xx} fc = l { CTl }fe = l 

where — \/(k&T), fee is Boltzmann's constant, T is 
the absolute temperature, the symbol denotes a 

summation over spin states of all the Ising spins and the 
symbol Tr^ stands for a trace over spin degrees of free- 
dom of the fcth Heisenberg trimer. After tracing out the 



spin degrees of freedom of the fcth Heisenberg trimer, the 
cluster partition function Zk will solely depend on the 
three nodal Ising spins Uf-i, (Tk2, and Cfc3 attached to the 
fcth Heisenberg trimer. Moreover, one may take advan- 
tage of the generalized star-triangle transformatio n 
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which can be used in order to replace this effective Boltz- 
mann's factor by the equivalent expression depending 
on the three nodal Ising spins only. The above map- 
ping transformation generally represents a set of eight 
algebraic equations, which can be obtained from Eq. (|5j) 
by substituting all available spin configurations of three 
Ising spins involved therein. In an absence of the exter- 
nal magnetic field, one actually gets just two independent 
equations; the first equation for two uniform configura- 
tions with the three equally aligned Ising spins 

Vi = Z fc (±l/2,±l/2,±l/2) 

= 2exp ( ^/3J H J cosh \%&J\ 
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and the second equation for six non-uniform configura- 
tions with one Ising spin pointing in opposite with respect 
to the other two 

F 2 = Z fe (±l/2,±l/2,Tl/2) = Z fc (±l/2,Tl/2,±l/2) 
= Z fe ( T l/2,±l/2,±l/2) 
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For abbreviation, we have introduced here two auxiliary 
functions Q defined as follows 
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It is noteworthy that two algebraic equations ©-© un- 
ambiguously determine yet unspecified mapping param- 
eters A and J c g 
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which ensure a general validity of the star-triangle trans- 
formation (JSJ on assumption that the mapping parame- 
ters A and J c g are given by Eqs. 

Substituting the generalized star-triangle transforma- 
tion (|S|) into the relation ^ one easily gets a rigorous 
mapping correspondence between the partition function 
Z of the spin- 1/2 Ising- Heisenberg model on the TIT lat- 
tice and respectively, the partition function Z\m of the 
spin-1/2 Ising model on a simple triangular lattice 



Z(fi, J H , Ji,A) = A~' N Zj M (p,>yJ eS ) ) 



(10) 



which is unambiguously given by the Hamiltonian with 
the effective nearest-neighbor interaction 7 J c ff 



Him = — jJeS ^ ofcrj 
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Despite the different spatial symmetries, the Ising- 
Heisenberg model on the two distinct TIT lattices can 
be rigorously mapped to the same effective Ising model 



on a simple triangular lattice, since the generalized star- 
triangle transformation (|5|) is performed locally for the 
same smallest commuting spin cluster in both the mod- 
els (see Fig. (Uc)-(d)). As a matter of fact, the most 
crucial difference between both the mappings lies just 
in a strength of the effective nearest-neighbor coupling 
7 J e ff, which is twice as large for the TIT lattice shown in 
Fig. [Hb) compared to that displayed in Fig. [TJa). The 
other mapping parameter A represents just less impor- 
tant multiplicative factor in the mapping relation (|10[) 
between both partition functions. It is worthwhile to re- 
mark, moreover, that the exact result for the partition 
function of the spin-1/2 Ising model on the triangular 
lattice is well knowni£~— 
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and hence, the partition function of the spin-1/2 Ising- 
Heisenberg model on the TIT lattices can also be exactly 
calculated by inserting the exact result (fT2"|) for the par- 
tition function Z\m of the spin-1/2 Ising model on the 
triangular lattice with the effective nearest-neighbor in- 
teraction 7 J c ff into the rigorous mapping relation (|10[) . 

At this stage, other basic thermodynamic quantities 
can also be calculated in a rather straightforward way 
within the framework of this rigorous mapping method. 
Using the standard relations of thermodynamics and sta- 
tistical physics, the free energy F of the spin-1/2 Ising- 
Heisenberg model on the TIT lattices can be expressed 
through the free energy Fjm of the corresponding spin- 
1/2 Ising model on the triangular lattice 



F = F 1M (/3 n J cS ) - 7 A^ B Tln A 



(13) 



Similarly, the internal energy U of the spin-1/2 Ising- 
Heisenberg model on the TIT lattices can be connected 
to the internal energy Uim of the spin-1/2 Ising model 
on the triangular lattice with the temperature-dependent 
nearest-neighbor interaction 7 J c ff 

4 V V 2 ) Jcff ^ Vx V 2 ) K ' 

For completeness, let us merely quote that the explicit 
form of the internal energy Ujm of the spin-1/2 Ising 
model on the triangular lattice can be for instance found 
in Ref . 49] , while two new functions W\ and W% entering 
the established mapping relationship (fl"4"|) for the internal 
energy are explicitly given by the following formulas 
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Furthermore, one may also readily obtain exact results 
for the entropy S and specific heat C by making use of 
the standard thermodynamic relations 



resultant exact expression for the single-site magnetiza- 
tion of the Ising spins in the spin-1/2 Ising-Heisenberg 
model on the TIT lattice 



k B lnZ + ^, 
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after substituting the previously derived exact results for 
the partition function (fpQ|) and the internal energy (fT4|) 
into the relations (fl6|) . 

Next, let us calculate the spontaneous magnetization 
of the Ising and Heisenberg spins, respectively. The spon- 
taneous magnetization of the Ising spins can be simply 
calculated by exploiting the exact mapping theorems de- 
veloped by Barry and co-workers^— 
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according to which the canonical ensemble average (• • •) 
of any function j\ of the Ising spins in the spin-1/2 Ising- 
Heisenberg model on the TIT lattice directly equals the 
canonical ensemble average (• • -)im of the same function 
fi of the Ising spins in the equivalent spin-1/2 Ising 
model on the triangular lattice. This result is taken to 
mean that the spontaneous magnetization of the Ising 
spins mi in the spin-1/2 Ising-Heisenberg model on the 
TIT lattice can in turn be calculated from the single-site 
magnetization mim of the spin-1/2 Ising model on the 
triangular lattice^ 
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with the parameter z = exp(— j/3J c s/2). Substituting 
the explicit form @ of the effective nearest-neighbor in- 
teraction 7 J e ff into the formula (TT5|) one finally gets the 
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whereas 7 = 1 for the TIT lattice shown in Fig. HJa) and 
7 = 2 for the TIT lattice drawn in Fig. \T[b)- 

Owing to the commuting character of different cluster 
Hamiltonians, the ensemble average of any function f% 
of the Ising and Heisenberg spins belonging to the kth 
cluster Hamiltonian ^ can be calculated from the gen- 
eralized form of Callen-Suzuki spin identity^— 
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The exact spin identity (|20[) allows us to express, after 
straightforward but a little bit tedious calculation, the 
spontaneous magnetization of the Heisenberg spins mn 
in terms of the triplet correlation t\M = (&ki&k2&k3)iM 
between three Ising spins from the fcth cluster Hamilto- 
nian and the formerly derived the single-site magnetiza- 
tion miM of the Ising spins 
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and two newly defined functions Qi and Q2 
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To complete our exact calculation of the single-site mag- 
netization of the Heisenberg spins mu it is now sufficient 
to find the triplet correlation t\u = (^ki^^^ks)™ of the 
spin-1/2 Ising model on the triangular lattice, which has 
been rigorously calculated by Baxter and Choy^ 
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where y = exp(7/3J e g) and z = exp(— 7/3 J c g /2) as be- 
fore. Substituting the effective nearest-neighbor interac- 
tion 7 Jog into the parameters y and z entering the for- 
mula (|23|) one obtains the required closed-form expression 
for the triplet correlation of the Ising spins 
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Last but not least, let us make a few comments on a 
critical behavior of the spin-1/2 Ising- Heisenberg model 
on the TIT lattices. It should be noted here that the 
critical behavior is always accompanied with a non- 
analytic behavior of the partition function or some of its 
higher-order temperature derivative. However, it may 
be easily understood from Eqs. ([5])-© that the map- 
ping parameter A as well as any temperature derivative 
of it is a smooth continuous function and therefore, it 
must not cause a critical behavior of the spin-1/2 Ising- 
Heisenberg model. Bearing this in mind, the spin-1/2 
Ising- Heisenberg model on the TIT lattices becomes crit- 
ical if and only if the partition function of corresponding 
spin-1/2 Ising model on the triangular lattice becomes 
critical as well. Accordingly, the critical condition al- 
locating critical points of the spin-1/2 Ising- Heisenberg 
model on the TIT lattices can readily be obtained from 
a comparison of the effective coupling of the correspond- 
ing spin-1/2 Ising model on the triangular lattice with its 
critical value 
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Here, it is emphasized that the inverse critical temper- 
ature (3 C — l/(fcB?c) enters the parameters V\ and Vi 



given by Eqs. (|6|) and ([7]) instead of f3. It is quite 
evident from the critical condition (|23|) that the spin- 
1/2 Ising- Heisenberg model on the TIT lattice is spon- 
taneously long-range ordered whenever the effective cou- 
pling is greater than the critical value 7/3 J c ff > In 3, oth- 
erwise it becomes disordered for 7/3J G ff < In 3. Thus, an 
intersection of the effective coupling "f(3 J e ff with its crit- 
ical value provides a feasible criterion for obtaining the 
critical points of the spin-1/2 Ising- Heisenberg model on 
both TIT lattices, which will be thoroughly discussed in 
Section [ETTB. 



III. RESULTS AND DISCUSSION 

In this section, let us proceed to a discussion of the 
most interesting results obtained for the ground state 
and finite-temperature properties of the spin-1/2 Ising- 
Heisenberg model on two considered TIT lattices. First, 
it is worth mentioning that all the results derived in 
the foregoing section are valid regardless of whether the 
interaction constants are assumed to be ferromagnetic 
or antiferromagnetic. While the change in character of 
the Heisenberg interaction has a profound effect upon 
the magnetic behavior of the spin-1/2 Ising- Heisenberg 
model on the TIT lattices, the respective change in the 
Ising interaction J\ — > — J\ merely causes a rather trivial 
spin reversal of — > —erf of all the Ising spins. There- 
fore, we will henceforth restrict our attention only to the 
particular case of the model with the ferromagnetic Ising 
interaction Jj > 0, while the respective behavior of the 
model with the antiferromagnetic Ising interaction J\ < 
can simply be deduced from the results to be presented. 
Let us finally mention that the Ising interaction will be 
subsequently used as the energy unit when defining two 
dimensionless parameters: the dimensionless tempera- 
ture kftT/Ji and a relative strength of the Heisenberg 
interaction with respect to the Ising one Ju/ Ji- 



A. Ground state 

Let us examine first the ground state of the spin- 
1/2 Ising- Heisenberg model on two geometrically related 
TIT lattices under consideration. Because of the com- 
muting character of different cluster Hamiltonians, the 
overall ground-state spin arrangement is simply given by 
a tensor product over the lowest-energy eigenstates of 
the cluster Hamiltonian If the lowest-energy eigen- 
value of the cluster Hamiltonian ^ is degenerate, how- 
ever, there does not exist a unique ground state but the 
highly degenerate ground-state manifold whose degener- 
acy might represent rather complex counting problem. 
Typical ground-state phase diagrams in the A — Jn/Ji 
plane including all possible ground states are depicted 
in Fig. [2Ia)-(b) for both investigated TIT lattices. Al- 
though the ground-state boundaries are identical for both 
the investigated models, there is a fundamental difference 
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in the character of spin arrangements of some ground- 
state phases. As one can see from Fig. [Ha), the Ising- 
Heisenberg model on the first TIT lattice from Fig. Ufa) 
displays two spontaneously long-range ordered and one 
disordered ground state, more specifically, the classical 
ferromagnetic phase (CFP), the quantum ferromagnetic 
phase (QFP), and the disordered quantum paramagnetic 
phase (QPP). All three phases are mutually separated by 
two lines of discontinuous (first-order) phase transitions 
given by the conditions 



CFP-QFP: 
CFP-QPP(ODP): 



Jh 
Ji 
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(for A > 1)(26) 
(27) 



which were obtained by comparing the energy of the 
lowest-energy eigenstates constituting the respective 
ground states. As could be expected, CFP with a perfect 
alignment of all the Ising as well as Heisenberg spins 
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dominates in the prevailing region of the parameter space 
with the ferromagnetic Heisenberg interaction Jh > 0. 
If the Heisenberg intra-trimer interaction has a strong 
easy-plane anisotropy, however, one may also encounter 
a more striking spontaneous long-range order 



TV N-y 

iQFP>=nit> CTf n^(i^) 

i=l fe=l* 



-int>+im»5f, 



£29) 



in which a perfect alignment of all the Ising spins 
is accompanied with a symmetric quantum superposi- 
tion of three up-up-down spin states of the Heisenberg 
trimers. Accordingly, the spontaneous magnetization of 
the Heisenberg spins undergoes a quantum reduction of 
the magnetization to one third of the saturation magne- 
tization and hence, we will refer to QFP as to the quan- 
tum ferromagnetic phase in view of a parallel orienta- 
tion of both sublattice magnetizations of the Ising and 
Heisenberg spins, respectively. An origin of this uncon- 
ventional spontaneous long-range ordering can be related 
to a competition between two different but ferromagnetic 
interactions, namely, the easy-axis Ising and easy-plane 
Heisenberg (A > 1) interaction. 

Last but not least, the parameter region with the an- 
tiferromagnetic Heisenberg interaction is dominated by 
QPP without any spontaneous long-range order. The 
disordered phase occurs on behalf of a spin frustration 
of the Heisenberg trimers, which are incapable to simul- 
taneously satisfy specific requirements of the antifcrro- 
magnetic intra-trimer interaction between the Heisen- 
berg spins. If the antiferromagnetic intra-trimer inter- 
action is stronger than the threshold value (27), then, 
the Heisenberg trimers have a tendency to form one sin- 
glet dimer and one unbound spin. The Ising spin coupled 
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FIG. 2: (Color online) Ground-state phase-diagrams of the 
spin- 1/2 Ising-Heisenberg model on the two TIT lattices 
shown in Fig. \V[ a) and (b). For a detailed description of 
the phases see the text. 



to two Heisenberg spins creating singlet pair will be con- 
sequently decoupled from the rest of the respective spin 
cluster and in this way, the spin frustration of the Heisen- 
berg trimers is also mediated to the Ising spins. Despite 
the spin frustration, the Ising spins still have a tendency 
of showing ferromagnetic short-range order as evidenced 
by the zero-temperature limit of the effective coupling 
lim/r^o 7/3 Jeff = ln2 (7 = 1) in the corresponding spin- 
1/2 Ising model on the triangular lattice implying that 
(cr?(r|) n . n , ~ 0.06564. The ferromagnetic short-range or- 
der of the Ising spins is related to tenfold degeneracy of 
the lowest-energy eigenstate of the cluster Hamiltonian 
(3), whereas four out of ten lowest-energy eigenstates of 
the Heisenberg trimer occur as long as the three enclosing 
Ising spins are equally aligned and another six provided 
that one out of three Ising spins points in opposite with 
respect to the other two (see Fig. [3] for pictorial repre- 
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If the three enclosing Ising spins are equally aligned (say 
all three Ising spins are pointing up), one actually finds 
just the two linearly independent lowest-energy eigen- 
states composed of one singlet dimer and one Heisenberg 
spin polarized by the two surrounding Ising spins into 
their direction (see the first four eigenstates in Eq. (|3"0")) 
schematically drawn in the first two rows of Fig. [3]) . On 
the other hand, there exists just one lowest-energy eigen- 
state of the Heisenberg trimer being composed of the 
singlet pair and one polarized Heisenberg spin on as- 
sumption that one Ising spin is pointing in opposite with 
respect to the other two. Under this condition, the po- 
larized Heisenberg spin is always connected to the two 
equally aligned Ising spins and the singlet pair decou- 
ples the third unequally oriented Ising spin from the re- 
maining spins of the relevant spin cluster (see the last 
six eigenstates in Eq. (|30)) schematically illustrated in 
the last two rows of Fig. [3]). Altogether, the ground- 
state manifold inherent to the QPP will be characterized 
by a rather huge but non-trivial macroscopic degeneracy, 
which will be explored in a more detail in Section IIIID 
by the analysis of entropy. 

The other disordered state DP emerging for the partic- 
ular limiting case of A = (i.e. the Ising limit) deserves 
a special mention, since there appears a dramatic change 
in the character of the lowest-energy eigenstates that con- 
stitute the ground-state manifold due to a complete lack 
of local quantum fluctuations. It can be easily proved 
that there exist just three different lowest-energy spin 
configurations on the smaller inner triangles for any given 
configuration of the three enclosing Ising spins (see Fig. fj] 
for the pictorial representation of the lowest-energy con- 
figurations) . Accordingly, the overall degeneracy of the 
ground-state manifold becomes for the DP phase a quite 
simple counting problem ft — 2 N 3~' N = (2 x 3 7 ) Ar , be- 
cause the Ising spins residing on the sites of the larger tri- 
angles contribute to the overall degeneracy by the factor 
2 N and the additional factor 3 7JV is related to the three- 
fold degeneracy of the lowest-energy spin configurations 
on the smaller inner triangles. Obviously, the similar 
procedure cannot be used for finding out the degener- 
acy of the ground-state manifold for the other disordered 
QPP of the more general Ising-Heisenberg model, since 




FIG. 3: The lowest-energy eigenstates of the cluster Hamilto- 
nian ((3]), which constitute the ground-state manifold for the 
QPP and ODP phases. The spin polarization is given by an 
arrow and an oval denotes a singlet-dimer state. 



there emerge either two or one lowest-energy eigenstate 
depending on whether the three enclosing Ising spins are 
equally aligned or not. 

As far as the ground state of the spin- 1/2 Ising- 
Heisenberg model on the other TIT lattice (Fig. Efb)) 
is concerned, the ground-state boundaries and the re- 
spective spin arrangements of two previously reported 
spontaneously long-range ordered phases CFP and QFP 
remain unchanged. However, it surprisingly turns out 
that the Ising-Heisenberg model on the second TIT lat- 
tice displays even in the highly frustrated region a quite 
peculiar spontaneous long-range order instead of the dis- 
ordered QPP for arbitrary but non-zero anisotropy pa- 
rameter A^O (note that the character of the other disor- 
dered state DP to emerge for A = remains unchanged). 
The existence of this unconventional spontaneous order- 
ing could be ascribed to the quantum order- from- disorder 
effect, which partially lifts the ground-state degeneracy 
of QPP due to a strengthening of local quantum fluc- 
tuations invoked by a higher number of the Heisenberg 
trimers in the second TIT lattice shown in Fig.QJb). The 
quantum order-from-disorder effect generally acts against 
the spin frustration and it may thus cause an appearance 
of the unusual partially ordered and partially disordered 
phase (ODP) with a spontaneously broken Z2 symme- 
try. Note furthermore that typical spin arrangements to 
emerge in ODP are completely the same as specified by 
Eq. (|3TJ)) for QPP and schematically illustrated in Fig. [3] 
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(a) 



FIG. 4: The lowest-energy spin configurations, which consti- 
tute the ground-state manifold for the DP phase to emerge 
in the Ising limit A = 0. Another twelve lowest-energy spin 
configurations can be obtained from the spin configurations 
drawn in the last two lines by rotating them around the three- 
fold C3 symmetry axis. 



The remarkable spontaneous order of the ODP will be 
convincingly evidenced in Sections IIIIB and MIC by a 
more comprehensive analysis of the critical behavior and 
the order parameter. However, the spontaneous order 
of ODP could also be evidenced by the effective cou- 
pling of the corresponding spin- 1/2 Ising model on the 
triangular lattice tending towards the low-temperature 
asymptotic limit lim/r^o 7/3J e ff = In 4 for 7 = 2, which 
is evidently above the critical value "ff3 c J c R = In 3 and 
is thus consistent with the spontaneously long-range or- 
dered ground state. It should be pointed out, moreover, 
that the finite zero-temperature limit of the effective cou- 
pling also serves in evidence of an imperfect spontaneous 
long-range ordering to emerge in ODP, which is also 
confirmed by non-zero albeit not fully saturated zero- 
temperature values of the spontaneous magnetizations of 
the Ising and Heisenberg spins m\ = |(yj§) 5 — 0-47482 

and m H = - §)(if§)* - 0.16408 obtained from 

Eqs. (|T9| and (|21"j) . respectively. Altogether, it could 
be concluded that the local quantum fluctuations are re- 
sponsible in ODP not only for the quantum reduction of 
the magnetic moment of the Heisenberg spins, but more 
strikingly, they also indirectly cause a small but non-zero 
(cca. 5%) quantum reduction of the magnetic moment 
of the classical Ising spins as well. 
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FIG. 5: (Color online) Temperature dependences of the effec- 
tive nearest-neighbor interaction 7/3J e g of the corresponding 
spin- 1/2 Ising model on the triangular lattice for the specific 
choice of the anisotropy parameter A = and various values 
of the interaction ratio Jh/Ji- Broken lines show the critical 
value of the effective coupling above (below) which the spin 
system becomes spontaneously ordered (disordered). 



B. Critical behavior 

Before proceeding to a detailed discussion of the crit- 
ical behavior, let us explore in detail temperature de- 
pendences of the effective interaction 7^ J e ff of the corre- 
sponding spin- 1/2 Ising model on the triangular lattice as 
depicted in Figs. [5] and [5] for two TIT lattices at various 
values of the interaction ratio Ju/Ji and two different 
values of the exchange anisotropy A = and 1, respec- 
tively. It can be clearly seen from Fig. [S] that the effec- 
tive coupling in the Ising limit A = either diverges for 
Jk/Ji > —1 or it tends towards zero for Ju/Ji < — 1 as 
temperature goes to zero. In the latter case, the effective 
coupling is always below its critical value "ff3 c J c ff = In 3 
and this result proves an existence of the DP in the Ising 
limit A = of both the investigated TIT models for 
any temperature if Ju/Ji < — 1- Contrary to this, the 
effective coupling 7/3J e ff either diverges or asymptoti- 
cally reaches some finite value for arbitrary but non-zero 
anisotropy parameter A ^ as shown in Fig. [S] for one 
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FIG. 6: (Color online) Temperature dependences of the effec- 
tive nearest-neighbor interaction 7/3J e s of the corresponding 
spin- 1/2 Ising model on the triangular lattice for the specific 
choice of the anisotropy parameter A = 1 and various values 
of the interaction ratio Jh/Ji- Broken lines show the critical 
value of the effective coupling above (below) which the spin 
system becomes spontaneously ordered (disordered). 



particular choice of A. While the effective coupling /3 J e g 
of the first TIT lattice always remains below its critical 
value in the highly frustrated regime Jh/Ji < —2/(2 + A) 
(Fig. [H[a)), the twice as large effective coupling 2f3J D g 
of the second TIT lattice is strong enough in order to 
induce the spontaneous long-range order at sufficiently 
low temperatures even in the highly frustrated region 
Jh/Ji < -2/(2 + A) (Fig. Mb)). This result thus pro- 
vides an independent confirmation of the quantum order- 
from-disorder effect, which arises from the local quan- 
tum fluctuations governing the magnetic behavior of the 
Heisenberg trimers. Finally, it is worthy to recall that 
an intersection of the effective coupling 7/3J e ff with the 
relevant critical point of the spin-1/2 Ising model on the 
triangular lattice can be straightforwardly used in order 
to locate critical points of the spin-1/2 Ising- Heisenberg 
model on the TIT lattices, which in fact represents the 
numerical solution of the critical condition (1231) . 

The critical temperature of two investigated TIT mod- 
els is plotted in Fig. [7] against a relative strength between 
the Heisenberg and Ising interaction for several values of 
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FIG. 7: (Color online) Critical temperature feTc/Ji of the 
spin-1/2 Ising-Heisenberg model on two TIT lattices from Fig. 
QJa)-(b) as a function of the interaction ratio Jh/Ji for several 
values of the exchange anisotropy A. 



the anisotropy parameter A. If one considers the partic- 
ular case of the easy-axis exchange anisotropy A < 1, the 
critical temperature of the first lattice model (Fig. E[a)) 
monotonically decreases with decreasing the ratio Jh/Ji 
until it completely vanishes at the ground-state boundary 
(f27|) between CFP and QPP. The highest critical temper- 
ature of CFP 



lim 



Jl 



111 



(2 + V3) 



0.75933 



(31) 



can be accordingly acquired in the asymptotic limit 
Jh/Ji - ► 00. On the other hand, the critical temper- 
ature exhibits a more striking non-monotonous depen- 
dence when assuming the easy-plane exchange anisotropy 
A > 1, because the sufficiently strong ferromagnetic 
Heisenberg interaction then favors a presence of uncon- 
ventional quantum ferromagnetic ordering QFP before 
the classical CFP one. The critical temperature there- 
fore tends just to one third of the previously reported 
asymptotic value 



hm 



Jl 



3 In 



(2 + Vtj 



0.25311, 



(32) 
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since the magnetic moment of all the Heisenberg spins in 
QFP is reduced by local quantum fluctuations to one 
third of their magnetic moment in CFP. Besides, the 
most conspicuous increase of the critical temperature 
upon lowering the interaction ratio Jh/Ji can be detected 
in a close vicinity of the phase boundary ([26j) between 
QFP and CFP, which is attributable to the increase of 
the magnetic moment of the Heisenberg spins when pass- 
ing from QFP towards CFP. As far as the critical behav- 
ior of the second TIT lattice (Fig. [7Jb)) is concerned, 
the critical temperature generally exhibits qualitatively 
the same dependences on assumption that the Heisen- 
berg interaction is ferromagnetic. As a matter of fact, 
the critical temperature is then shifted towards slightly 
higher values due to a higher connectivity of the Ising 
spins in the other TIT lattice and the asymptotic values 
of the critical temperature in Jh/ Ji oo limit read 



lim 



feB^c 

, Ji 



In 



1 + V3 + V2\/3 



for CFP if A < 1 and 

&bJ c 



lim 



Ji 



3 In 



^ 1 + V3 + V273 ^ 



1.20273 (33) 



0.40091(34) 



for QFP if A > 1. The most fundamental difference 
in the critical frontiers of two investigated TIT mod- 
els can be thus found in the parameter region with a 
strong (negative) antiferromagnetic Heisenberg interac- 
tion, where the critical lines of the latter model do not 
vanish but they tend towards some non-zero value af- 
ter passing through a global minimum in a proximity of 
the phase transition ([27)) between CFP and ODP phases. 
It could be concluded that the presented exact results 
for the critical boundaries provide another independent 
confirmation for an existence of the striking spontaneous 
long-range order ODP, which emerges in the latter TIT 
model in spite of the high spin frustration. 



C. Spontaneous magnetization 

Now, let us turn our attention to a discussion of typical 
temperature dependences of the spontaneous magnetiza- 
tion, which will bring insight into the main differences in 
the thermal behavior of three spontaneously long-range 
ordered phases CFP, QFP, and ODP. For this purpose, 
the temperature dependences of both spontaneous sub- 
lattice magnetizations are plotted in Figs. [5] and [§] for 
two investigated TIT lattices, two different values of the 
exchange anisotropy A and several values of the inter- 
action ratio Jh/Ji- First, let us take a closer look at 
thermal variations of the spontaneous magnetizations of 
the Ising and Heisenberg spins in two TIT lattices with 
the isotropic Heisenberg intra-trimer interaction. It is 
quite obvious from Fig. HJa) that the spontaneous mag- 
netizations of the Ising and Heisenberg spins for the first 



TIT lattice start from their saturated values, which bear 
evidence of CFP unless the disordered QPP becomes the 
ground state for Jh/Ji < —2/3. In addition, it is also 
quite clear from Fig. Ufa) that the spontaneous mag- 
netizations of the Ising and Heisenberg spins exhibit a 
similar temperature-induced decline, the magnetization 
of the Heisenberg spins actually shows only a slightly 
greater temperature-induced downturn than the sponta- 
neous magnetization of the Ising spins even though both 
sublattice magnetizations tend to zero with the same 
critical exponent f3 c = 1/8 from the standard univer- 
sality class of the two-dimensional Ising model. It should 
be nevertheless noted here that the same general trends 
can also be detected in the respective temperature de- 
pendences of the spontaneous magnetizations of the sec- 
ond TIT lattice whenever J H /Ji > -2/3 (see Fig.^b)). 
However, the spontaneous magnetizations of the Ising 
and Heisenberg spins in the latter TIT model start from 
non-zero but not fully saturated values mi ~ 0.47482 
and 771h — 0.16408 even if Jh/Ji < —2/3, which bear ev- 
idence of the unconventional spontaneous ordering ODP 
basically affected by the local quantum fluctuations. The 
spontaneous magnetization of the Heisenberg spins re- 
veals the transition between CFP and ODP through a 
vigorous temperature-induced increase of this sublattice 
magnetization as long as the interaction ratio Jh/Ji is se- 
lected sufficiently close to but slightly below the ground- 
state boundary ([27]) (see the curves for Jh/Ji = —0.7 
and -0.75 in Fig.^b)). 

To shed light on another unusual spontaneous long- 
range ordered phase QFP, Fig. [HI displays temperature 
variations of the spontaneous magnetization of two inves- 
tigated TIT lattices for the ferromagnetic Heisenberg in- 
teraction with the easy-plane exchange anisotropy A = 2. 
Under this condition, the spontaneous magnetizations of 
the Ising and Heisenberg spins apparently exhibit accord- 
ing to Fig. [9] qualitatively the same temperature depen- 
dences for both investigated TIT lattices. If the rela- 
tive strength between the Heisenberg and Ising interac- 
tion exceeds the boundary value ([26]). one actually finds 
that the initial value of the spontaneous magnetization 
of the Heisenberg spins is reduced by the local quantum 
fluctuations to one third of its saturation magnetization 
while the initial value of the spontaneous magnetization 
of the Ising spins is still fully saturated. This observation 
is fully consistent with the ground-state spin arrange- 
ment that was attributed to QFP. Beside this, thermal 
excitations in QFP may give rise to an interesting in- 
crease in the spontaneous magnetization of the Heisen- 
berg spins, whereas the observed temperature-induced 
increase in mn is the greater, the closer the interaction 
ratio Jh/Ji is selected to the ground-state boundary (j26[) . 
Note furthermore that one may also get the special value 
of the spontaneous magnetization of the Heisenberg spins 
toh = 1/3 due to the coexistence of the CFP and QFP if 
the interaction ratio is selected exactly at their ground- 
state boundary l|26p. i.e. Jh/Ji = 1 for the particular 
case of A = 2. 
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FIG. 8: (Color online) Temperature dependences of the spon- 
taneous magnetization mi of the Ising spins (broken lines) and 
the spontaneous magnetization iiih of the Heisenberg spins 
(solid lines) for the fixed value of the anisotropy parameter 
A = 1 and several values of the interaction ratio Jh/Ji. 
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FIG. 9: (Color online) Temperature dependences of the spon- 
taneous magnetization mi of the Ising spins (broken lines) and 
the spontaneous magnetization mu of the Heisenberg spins 
(solid lines) for the fixed value of the anisotropy parameter 
A = 2 and several values of the interaction ratio Jh/Ji- 



D. Entropy 

Next, let us turn to a detailed analysis of the entropy, 
which enables a deeper insight into the degree of random- 
ness in the disordered states. Figs . [TUl and ITT1 depict tem- 
perature variations of the entropy per one spin S'/A^t^b 
(iVr = N(l + 37) is the total number of all spins) for two 
different values of the exchange anisotropy and several 
values of the interaction ratio Jh / J\ ■ It is quite evident 
from Fig. [TU] that the highly macroscopically degenerate 
ground state DP develops in the Ising limit A = of 
both investigated lattice models due to the spin frustra- 
tion, which originates from the sufficiently strong antifer- 
romagnetic intra-trimer interaction Jh/Ji < — 1. Con- 
sequently, the DP manifests itself through a relatively 
large residual entropy S/NtI^b = jln.6 ~ 0.44794 and 
I In 18 ~ 0.41291, which is in agreement with the calcu- 
lated degeneracy of the ground-state manifold inherent 
to the DP in the limiting Ising case (see Section IIIIR ) . 
On the other hand, arbitrary but non-zero anisotropy 
parameter A 7^ is responsible for an onset of the lo- 
cal quantum fluctuations, which lift the macroscopic de- 



generacy approximately by 50% in the highly frustrated 
regime of both TIT lattices (see Fig. Hi] for A = 1). 
While the residual entropy of the former TIT lattice 
SyiVrfcB — 0.23133 is still high enough to preserve the 
disordered nature of QPP, the slightly lower residual en- 
tropy of the latter TIT lattice S/N T k B ^ 0.20075 sup- 
ports an existence of partially ordered and partially dis- 
ordered state ODP. In spite of the spontaneously broken 
symmetry, the ODP still retains a relatively high degree 
of randomness and is quite reminiscent of the fully dis- 
ordered QPP. 



E. Specific heat 

Finally, let us discuss the main features of tempera- 
ture dependences of the zero-field specific heat, which 
are displayed in a semilogarithmic scale in Figs. [12] and 
Q2] for two investigated TIT lattices by selecting two dif- 
ferent values of the exchange anisotropy A and various 
values of the interaction ratio Jh/Ji- Typical thermal 
variations of the specific heat for the particular case of 
the isotropic Heisenberg intra-trimer interaction are il- 
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FIG. 10: (Color online) Semilogarithmic plot for tempera- 
ture dependences of the entropy per one spin S/NtIcb for 
two considered TIT lattices, the special value of the exchange 
anisotropy A = and several values of the interaction ratio 
Jh/Ji. 
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FIG. 11: (Color online) Semilogarithmic plot for tempera- 
ture dependences of the entropy per one spin S/Nt1cb for 
two considered TIT lattices, the special value of the exchange 
anisotropy A = 1 and several values of the interaction ratio 
Jh/Ji. 



lustrated in Fig.Q^l As one can see from Fig. [T27a). the 
Ising-Heisenberg model on the first TIT lattice exhibits a 
logarithmic singularity from the standard Ising universal- 
ity class when considering the ferromagnetic Heisenberg 
interaction Jh > 0. It is quite evident from this fig- 
ure, moreover, that a gradual strengthening of the anti- 
ferromagnetic intra-trimer interaction causes an appear- 
ance of a shoulder superimposed on the low-temperature 
tail of the specific heat divergence (see the curve for 
Jh/Ji = —0.5). If the interaction ratio Jh/Ji is se- 
lected sufficiently close but slightly above the ground- 
state boundary (f2"Tf between CFP and QPP, the marked 
round maximum at relatively low temperatures is sub- 
sequently followed by the logarithmic divergence super- 
imposed on ascending part of another round Schottky- 
type maximum emerging at higher temperatures (see the 
curve for Jh/Ji = —0.65 in the inset of Fig. [TSTa)). In 
agreement with the ground-state and finite-tempcraturc 
phase diagrams displayed in Figs.[5Ja) and[7][a), the loga- 
rithmic singularity completely vanishes from the temper- 
ature dependence of the heat capacity with regard to the 
disordered character of QPP whenever the antiferromag- 
netic intra-trimer interaction exceeds the ground-state 
boundary (|27p between CFP and QPP (see the curves for 



Jh/Ji < —2/3). After passing through this ground-state 
boundary, the low-temperature round maximum rather 
steeply diminishes within QPP as the Heisenberg interac- 
tion further strengthens (see the curve for Jh/Ji = —0.68 
in the inset of Fig. HHa)), then it shows a single round 
maximum of a rather irregular shape for the Heisenberg 
interaction of a moderate strength (e.g. for Jh/Ji = — 1) 
and finally, the specific heat displays two more or less 
separated round maxima for a relatively strong Heisen- 
berg interaction (e.g. for Jh/Ji = —2). As far as the 
specific heat of the second TIT lattice model is con- 
cerned (Fig. [HJb)), it generally shows similar thermal 
variations except that one still encounters a marked log- 
arithmic singularity inside the highly frustrated region 
(Jh/Ji < —2/3 for A = 1), which is occupied by the un- 
conventional partially ordered and partially disordered 
state ODP rather than the fully disordered QPP. 

Last, let us briefly comment on a typical thermal be- 
havior of the heat capacity in a close vicinity of the an- 
other possible phase transition between CFP and QFP. 
For this purpose, Fig. [TJJ] illustrates typical temperature 
dependences of the specific heat for the particular case 
of the ferromagnetic Heisenberg interaction with one se- 
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FIG. 12: (Color online) Semilogarithmic plot for temperature 
dependences of the specific heat per one spin C/NtIzb for 
two considered TIT lattices, the special value of the exchange 
anisotropy A = 1 and several values of the interaction ratio 
Ja/Ji. 
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FIG. 13: (Color online) Semilogarithmic plot for temperature 
dependences of the specific heat per one spin C/NtIcb for 
two considered TIT lattices, the special value of the exchange 
anisotropy A = 2 and several values of the interaction ratio 
Jh/Ji- 



IV. CONCLUDING REMARKS 



lected value of the easy-plane exchange anisotropy A = 2. 
According to these plots, the specific heat of both consid- 
ered TIT lattices exhibits qualitatively the same thermal 
variations with only a small shift of the logarithmic sin- 
gularity of the latter model towards higher temperatures. 
If the interaction parameters are tuned sufficiently close 
to the ground-state boundary between CFP and QFP 
(Jh/Ji = 1 for A = 2), then, one observes a develop- 
ment of the round Schottky-type maximum at relatively 
low temperatures reflecting spin excitations from CFP 
to QFP (for Jh/Ji < 1) or vice versa (for Jh/Ji > 1). 
The closer is the ratio Jh/Ji to the ground-state bound- 
ary (|26[) . the more pronounced the low-temperature 
Schottky-type maximum can be observed in the relevant 
thermal dependence. In addition to this rather robust 
low-temperature round maximum, there also may appear 
the second high-temperature round maximum upon fur- 
ther increase of a relative strength of the Heisenberg in- 
teraction as convincingly evidenced by temperature vari- 
ations of the specific heat displayed in the inset of Fig. Q2] 
for Jh/Ji = 1-3. 



In the present article, the spin- 1/2 Ising- Heisenberg 
model on two different but geometrically related TIT 
lattices has been exactly solved through the generalized 
star-triangle transformation establishing a rigorous map- 
ping equivalence with the corresponding spin- 1/2 Ising 
model on a triangular lattice. Within the framework 
of this rigorous mapping method, we have derived ex- 
act analytical results for several basic thermodynamic 
quantities such as the free and internal energy, sponta- 
neous magnetization, entropy, specific heat, and we have 
also constructed the ground-state and finite-temperature 
phase diagrams quite rigorously. It has been demon- 
strated that the spin-1/2 Ising-Heisenberg model on two 
TIT lattices exhibit a surprisingly rich magnetic behav- 
ior including several unconventional and yet undetected 
quantum phases without any classical counterpart. 

A mutual competition between two ferromagnetic in- 
teractions of basically different character (the easy-axis 
Ising and easy-plane Heisenberg interactions) is the main 
cause for the emergence of QFP in which a symmet- 
ric quantum superposition of three up-up-down states 
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of the Heisenberg trimers accompanies a perfect align- 
ment of all the Ising spins. Apart from this common 
feature of both studied lattice models, it has been evi- 
denced that two investigated TIT lattices display a very 
different spin ordering (QPP versus ODP) in the highly 
frustrated region owing to the crucial difference in a rel- 
ative strength of local quantum fluctuations. Among the 
most remarkable findings one should mention an exis- 
tence of ODP alone, which is in part spontaneously or- 
dered as evidenced by the singular behavior and power- 
law decay of several thermodynamic quantities in a close 
vicinity of the critical point, and in part disordered, as 
evidenced by the non-zero residual entropy. It actually 
turns out that the local quantum fluctuations are respon- 
sible in ODP not only for a quantum reduction of the 
spontaneous magnetization of the Heisenberg spins, but 
they also indirectly cause a quite peculiar quantum re- 
duction of the spontaneous magnetization of the other- 
wise classical Ising spins. The present work thus provides 
the first example of exactly solved frustrated spin model 
with a mutual coexistence of the imperfect spontaneous 
order and the partial disorder, which still exhibits a non- 
trivial criticality at finite temperatures. To the best of 
our knowledge, this finding represents a quite novel phe- 
nomenon not reported in the literature so far. 

Finally, let us make a few remarks on the relevant mag- 
netic behavior of the spin- 1/2 Ising- Heisenberg model 
on the TIT lattices in the non-zero external magnetic 
field. While the star-triangle transformation can be 
rather straightforwardly adapted in order to account for 
the external magnetic field^ 5 - it will consequently estab- 
lish a rigorous mapping correspondence with the equiva- 
lent spin- 1/2 Ising model on a triangular lattice including 
the effective nearest-neighbor interaction along its edges, 
the effective triplet interaction within its faces and the 



effective field acting on its sites. Although the latter 
classical Ising model is in general not exactly tractable 
at finite temperatures, one still may take advantage of 
some precise numerical method (e.g. Monte Carlo sim- 
ulation) to examine the equivalent purely classical Ising 
model in an attempt to extract the magnetic behavior 
of the Ising-Heisenberg model from this exact mapping 
correspondence. In our follow-up paper we have there- 
fore restricted our attention to a rigorous analysis of 
the zero-temperature magnetization process of the Ising- 
Heisenberg model on the two considered TIT lattices, 
which is also compared with the relevant magnetization 
process of the analogous but purely quantum Heisen- 
berg model treated by means of the exact numerical 
diagonalization.— Our recent results show a remarkable 
coincidence between the zero-temperature magnetization 
processes of the Ising-Heisenberg and Heisenberg models 
displaying the same sequence of intermediate magnetiza- 
tion plateaus, which proves the usefulness of the simpli- 
fied Ising-Heisenberg model in providing guidance on a 
more complex magnetic behavior of the quantum Heisen- 
berg model. 59 ' 60 
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